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Abstract 

We derive and analyze the 'SHE' (Spherical Harmonics Expansion) type system of equations coupled in 
energy. We also show that diffusive behavior occurs on long time and distance scales and we determine the 
diffusion tensor. The analysis is based on the governing kinetic equations arising in electron transport in 
semiconductors. 
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1 Introduction 

This paper presents an extension of a previous work by P. Degond and S. Mancini in which a diffusion 
model describing the evolution of an electron gas confined between two parallel planes by a strong magnetic 
field is derived. Electrons colliding against the plane are supposed to be reemitted following a combination 
of diffusive and specular elastic laws. When the distance between these planes goes to zero, the distribution 
function (solution of a Vlasov-Poisson equations) is proved to converge to the solution of the macroscopic model, 
the so-called SHE model which is a diffusion model for the energy distribution function. This situation arises in a 
certain kind of ion propellers for satellites. SHE is an acronym for Spherical Harmonics Expansion coming from 
its earlier derivation by physicists |16l . These models appeared to be a good compromise between a very accurate 
description of physical phenomena by kinetic models and less numerically expensive macroscopic models such as 
Drift-Diffusion or Energy- Transport models. Whereas kinetic equations deal with functions of seven variables 
(six variables in the phase and time) and classical macroscopic models deal with functions of four variables (three 
in the position, space and time), SHE models introduce an intermediate one dimensional variable replacing the 
velocity and which appears as the kinetic energy associated to the velocity of a particle. This makes of SHE 
models quite reliable models at a rather low numerical cost. The SHE model in literature has been derived from 
the Boltzmann equation first by P. Dmitruk, A. Saul, and L. Reyna |16| . 

The mathematical theory of the diffusion approximation started with the seminal papers [HI, in the 
context of neutron transport after the formal theory had been set up by Hilbert, Chapman, Enskog and co- 
workers (sec for instance QJ). Their approaches were later extended by many authors. In this work, we are 
interested in a physical situation where the observation length scale is large compared to the mean free path 
while the observation time is large compared to the characteristic time evolution of the particle. 

The mathematical study of the vanishing mean free path limit and of diffusion approximation is by now a 
classical problem with applications in various fields of physics. We refer among others to 0, |24) and for recent 
|15j and for SHE models 0, |14| and the references therein for details as the physical background concerning 
these models. Concerning SHE models this study has been derived from the Boltzmann equation first by N. Ben 
Abdallah and P. Degond |7] under the assumption that the dominant scattering mechanisms are elastic collisions. 

*e-niail; dogbe@math.unicaen.fr 
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The present work has been inspired by 0, |14| . but differs from them in the nature of coUisions 
of electrons with the waU and the presence of a magnetic field. In |14| and |13) . elastically diffusive collisions 
at the plates is considered: the particles are reemitted with the same energy as their incident one, and with 
a random velocity direction. The force field is the gradient of a smooth potential function which is assumed 
given, independent of time and of the mean free path, and amounts to assume vary over the macroscopic scale 
only. As a consequence, the large time behavior of the distribution function is given by an energy distribution 
function depending on the parallel components relative to the boundary of the domain and satisfying a diffusion 
equation in both position and energy. Here, we treat the self-consistent problem, the electric potential satisfying 
the Poisson equation. Hence, in contrast of 1141 and 1131 . the presence of a quadratic term in Vlasov equation 
gives rise to some singular term which adds additional technical difficulties. The existence and regularity of 
the solution to the asymptotic model which is now constituted of the diffusion equation coupled to the Poisson 
equation, are not immediate, nevertheless verified. Our main contribution in this paper is to give a rigorous proof 
of this convergence. 

In p] the authors showed that collisions with a boundary can drive a system of neutral particles towards a 
diffusion regime and the diffusivity is infinite. It has been shown in 1151 that a logarithmic time re-scaling restores 
a finite diffusivity. Nevertheless the divergence which appears in [2| and waived in 1151 does not appear in our 
work because of the presence of a strong magnetic field directed transversally to the plates. Particle motion then 
does not occur along straight lines, but rather along the plates whose axes arc parallel to the magnetic field lines. 

We note also, that some results in the direction of this paper have been obtained, without the boundary 
conditions, and where collisions are taken into account through the non linear Pauli operator by T. Goudon et 
al in EHl. 



2 Setting up the problem and the main result 

Let us recall the physical background carefully. Our starting point is a Vlasov equation for the electron 
distribution. The microscopic model describes the evolution of particles limited between two parallel plates. In 
the region of this space, the particles are submitted to a force field. The physical space where the electrons 
evolve is £7 X R'^ where £7 = [0, 1] x R^. The position vector is denoted hy X = {x,y,z) and we split X into 
its perpendicular x e [0, 1] and parallel components ^ = {y, z) e relative to the boundary F. We denote by 

V = (vxjVyjV^) e R'^ the position and the velocity vectors of an electron between the planes. We decompose 

V = {vx,v), where is the velocity component parallel to the x-axis and v = {vy,Vz) £ R^ is the component 
parallel to the plates. The electrons are subject to a magnetic field transverse to the plates {B{(^), 0, 0) depending 
upon ^ and to a potential force field parallel to the plates, (0, Ey{^),E2{^)) depending only upon ^ and satisfying 
E_ = —\/^<l>, where V{ denotes the 2-dimensional gradient with respect to ^ and (f) is the potential, (cf. Fig.l). 
Therefore, they are supposed to move between the plates according to a collisional transport equation. 



E{y,z).e^ 







Fig. 1. Geometry of the model 



Here f{X,v,t)dXdv is understood to give the number of electrons that occupy any infinitesimal volume dXdv 
at the point (X, v) which is given as the solution of the boundary value problem: 



+ (li • + £ ■ V,) / + 1^ - (t; X B) • V„/ = 

where {vxB) = (0, v^B, —VyB), while the electrostatic potential <j) solves the Poisson equation 

f 1 



S = -V, 



fdxdv-C{^). 



(2.1) 



(2.2) 



We shall assume that C{^) is a given (nonnegative) function which is regular, say C°°(n). 
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Lot us specify the proper bouirdary conditions to be considered. We iirtroduce the set O = x R'^, aird 
consider its boundary F = 7 x R^, wliere 7 = {0, 1} x R^; we denote Rj. = {w £ R'^ : ±Vx > 0} and the 
following incoming and outgoing subsets of F (respectively representing incoming and outgoing particles to the 
domain Q,): 

T_ = ({0} X R''* X R^ U ({1} X R-' X R=L) 

where for instance, ({0} x R'^ x Rij_) represents electrons entering the region < x < 1, through the boundary 
plane x = Q. The boundary F_|_ is obtained by reversing the inequalities. We introduce the traces (or boundary 
values) of / on F x R^ according to: 

7(/) = /l{.=o.i} . 7*(/) = /l{.=o,i, ±..>o} : (2-3) 

7+(/) is the distributional function of the particles exiting the domain Q, at the boundary F, while 'y^{f) is that 
of incoming particles. We suppose, finally that it is a function of outgoing trace through an operator IC which 
expresses the interaction of the particles: 

7-(/)=^(7+(/)). (2.4) 
Introducing the operator B{^{f)) = 7^(/) — IC{^^{f)), we have B{^{f)) = 0. The boundary condition for (j) is: 

lim (t>{^,t) = 0, a.e. t > 0. (2.5) 

The homogeneous boundary condition for (j) means that the system of electrons is in equilibrium at infinity. We 
refer to for further properties of these boundary conditions and physical interpretations. In this work some 
important relations valid for these boundary conditions are derived. 

Our model is based on the assumption that the distance between the plates is small compared to the char- 
acteristic length of the horizontal motion. Since we are looking for a diffusion process, this suggests to rescale 
the longitudinal coordinate according to x' = ax, where a = 0(1) is the small parameter. But then, the number 
of collisions between a typical particle and the plates per unit of unsealed time is large: if we assume that the 
collisions between the particles and the plates are a purely isotropic process, there is no reason to expect that, 
at a large scale, the particles would follow a horizontal drift. Hence, in order to observe a horizontal motion at a 
large scale, it is logical to rescale the time variable as t' = a^t. After rescaling, (dropping the primes for the sake 
of clarity), the equation l|2.1|l - ^2.2|l is recast as: 

adtr + {v-yi + E'' • V„) /" + i (^v., ^-(vxB)- ,r = 0, (2.6) 

^- = -V50", -Air= [ [ rdxdv-C{Q. (2.7) 
It is clear that these scalings do not induce any modification in the boundary condition, namely, /" still satisfies: 

7"(r) = ^(7+(r)). (2.8) 

Finally, we prescribe an initial condition which is compatible with the expected asymptotic dynamics. This initial 
condition is 

f{x,v,t = o) = fi{x,v) y{x,v)enxR^. (2.9) 

Now, since we want to base our approximation of diffusion in position-energy space, it is relevant to interpret the 
velocity of particle in terms of its energy. For that purpose, let us introduce the spherical coordinates in velocity 
space: 

fLu = v/\v\ E S'^ (angular variable); 
e = |t;p/2, (energy variable); (2.10) 
V = \v\lu = spTtuj = N{t)uj. 

We identify v with the pair (e,ci;). Hence, for C"^ function e : R'^ ^ R, and any integrable function !^ e R'^ we 
have (coarea formula ^Tj): 

I (p{v)dv = [ [ ip{V2^uj){2e)^^'^d€dij (2.11) 
Jvev.^ Jo Js'-' 

(p{e,Lu)N{e)dedu} 
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where duj denotes the normalized Euclidean measure on S , so dv = r drdcu, with r = = v2e, and dr = 
In the coarea formula H2.1H . Af(e) represents the energy-density of states. 

The main contribution of this paper is to investigate the limit q ^ of 12. 61 - 12. 98 . We show that the limit 
/'^ as a — > is a function F{£^,e,t) of the longitudinal coordinate 5, of tlic energy e = |i'P/2 and of the time 
which obeys a diffusion equation in the position-energy space. In order to simplify the exposition, from now on, 
we will denote by (P) the problem JmJ, ^J^, (|2|HJ, EIJ. 
We are now ready to state our main result. 

Theorem 2.1 1. As a goes to zero, the solution /" to the problem (P) formally converges to an equilibrium 
state -F(5, |'(jp/2, f), solution to the following SHE model: 

BE { f) \ 

4^%/2i— + (s/i-E-j -1 = 0, (2.12) 

J(e,e,t) = -D(5,e) fvj-^l-V. (2-13) 



E = -V</., -A^= [ [ f{^,v,t)dxdv~C{0, (2.14) 
Jo 

E{^,e,t = 0) = Fi{^,e), (2.15) 

J(C,e = 0,t) = 0, (2.16) 

where Ej is a suitable initial condition, in the domain (^,e) G x (0, oo). The 'diffusivity tensor' H is 
given by: 



(2e)^/2 / / x{x,uj-^,e)®uj,dxdu, (2.17) 



I p 

where lo = {ujy,u!z), x = iXy^Xz), X®^ i/*^ tensor product {Xi'^j)i.je{y.z} o-'nd Xii^^i^'^i^-^) o- 
solution of the problem 

^Xi 

-'f.x-^ + (u X B) • V„Xi «" 6 

(2.18) 

7+(xO=^*(7"(xO) onV, 

satisfying i i xdxduj = 0, for all (^,£), where K.* is the operator adjoint of fC. 
Jo Js2 ~ ' ~ 

2. Under hypothesis 15. il (to be specified later on) and 15. gl there exists a solution {f",E^) in 
L°°{G,T;L'^{Q)) n 7/^(0, T; Lf„^(K|)) for any T e R+. f" converges to f° m the weak star topology of 
L°°{0,T-L'^{e)) for any T > 0, where f°{t,^,£) = E{t,^,£) and F{t,^,e) is the weak solution of the 
problem f^i - f^) . 

The result is very close to |T1]. Nevertheless, the proof will be different, though we shall use a lot of results 
developed in |14| . Existence results of Cauchy problem for Vlasov-Poisson system arc now well known (see |27) 'l 
in dimension three without external potential and the boundary value problem for Vlasov-Poisson was studied 
in 1221. The theory of global weak solutions is due to The detailed mathematical study of (P) in the case 
of a constant magnetic field and a given electrostatic potential, which varies only over the macroscopic scale will 
be found in |14|. We shall not dwell on the existence and uniqueness of a solution for the Cauchy problem (P) 
which can be done by means of Leray-Schauder's fixed point Theorem; we shall solely give a material which 
makes it possible to obtain this solution and focus on the establishment of the limit model. The investigation 
of this limit when a —f proceeds in two steps. The first one consists in showing that formally converges 
to a function of f(x,e,t) only, solution of diffusion 'SHE' problem. Then we derive the continuity and current 
equation 12.1211 , 12.131 . The second one corresponds to a rigorous convergence proof; we show that weak solutions 
of Vlasov-Poisson equations converge weakly (in an sense) towards weak solutions of the 'new' SHE model. To 
achieve these goals, two methods can be developed: the Hilbert expansion method [E| and the moment method 
[15| . We shall choose the latter because it involves more straightforward computations. 

The outline of the paper is as follows. Section is devoted to the preliminary materials regarding the 
functional setting of our problem and important properties of the boundary collision operator. We establish 
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some mathematical results on this operator. In particular, we prove that the flux of particles is conserved at the 
boundary and there is an entropy dissipation (that is, a Darrozcs-Guiraud-like inequality in gas surface interaction 
(121 . I^l). This will be enough to allows us to derive the formal asymptotic limit which is achieved in Section 
^ This asymptotic limit appears as a singular perturbation problem for Vlasov-Poissoir equations. We solve 
this problem and derive our SHE model as the corresponding limit equation. A necessary step for the rigorous 
derivation is to establish the existence of a solution for the kinetic problem; this will be done in Section|3 Finally 
Section El is dedicated to the rigorous proof of convergence itself. 

3 The boundary operator: assumptions and properties. 

In the sequel, we consider an expression of K, as follows 

K.^{v)= ( K{i,\v\^/2-uj' ^uj)(i){\v\u;')\uj'J,duj', Vt; £ R^ (3.1) 

i{c^'6S2, (|,'i>)er+} 

where u = {LUx,i^y,i^z) = = 1 is the decomposition of v into spherical coordinates, §^ = {oj e R'^, \uj\ = 

1} is the unit sphere. The operator IC maps the outgoing trace 7^(/) to the incoming one 7^(/) . The expression 
13.111 models an elastic bounce against the pianos with a random deflection of the velocity. Note that K is an 
integral kernel which describes the reflection law of the velocity direction. The quantity if (^, |f P/2; a;' — > Lj)|[i;j;|da; 
is the probability for a particle hittiirg the planes at point ^ with velocity v' = \v\lj' to be reflected with the same 
\v\ and velocity direction u in the solid angle duj. Bearing in mind that IC = IC(S,, \v\'^/2) operates on the angular 
variable uj only while ^ and \v\ are mere parameters. 

We give the adapted functional setting for the study of the collision operator. We set S±{x), a; = 0, 1 to be 
the following half-spheres: 

5+(0) =5_(1) = {a; e §^ o;^: < 0}, 5_(0) = 5+(l) = {l^ £ S^ o;^. > 0}. (3.2) 

We introduce the domain 5 = [0, 1] x with its associated incoming and outgoing boundaries defined by: 

5_ = (5_(0) X {0})U(5_(1) X {1}), 5+ = (5+(0)x{0})U(5+(l)x{l}). (3.3) 

We denote by L'^{S±) the space of square integrable functions on S± with respect to the measure |a;j,.|da;; by 
(/)5)e, (/, 3)r± the irmer products on 1/^(6) and on L^(r*) respectively defined by: 



(,Aff)e= / fgde, {f,g)r± = / fg\nAdT 

where d9 = dxd^^dv is the volume element in phase space, and dV = X]i:=o i d£,dv is the surface element and by 
(/i5)5i (/)5)s± the inner products on L'^{S) and L^(S±) defined analogously: 



(/.5)s = / / {fg){x,uj)dxduj, 

(/>5)5± = / / {fg){x,uj)\uj.x\dxdu;+ / / {fg){x,u;)\u;x\dxduj. 
JoJs±(o) JoJs±(i) 

and \ f\s, \.f\s± the associated norms. We shall denote also ©' = x R+* with dO' = dxde its volume element. 
V/ = /(a;), (/) will denote the angular average of / on the sphere and with respect to x i.e. 



if) = -r f{x,u))dxdu). 



We denote by = {/ e I/^(Rj_) : /(e,a;) is constant with respect to lu}. We define the operator as the 
orthogonal projectioir (for inner product (•, •)s^) of 1/^(5+) on the space C*, i.e. 

Q^f{x,Lu) = - f f{Lj)\Lj,\dw, LjeS±, xe{0,l}, (3.4) 

and the operator P*, as the orthogonal complement of : -P* = / — Q*. 

We shall list the required properties of the reflection operator K. They are summarized in the following 
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Hypothesis 3.1 (i) Flux conservation: 

I K{uj' uj)\LOj:\du) = 1, 

JS-(x) 



(3.5) 



for almost all {u),ijj') e S-{x) x S+{x), a; = 0, 1. 

(ii) Reciprocity principle: 

K{uj' uj) = K{-LO -Lo'), V(a;,a;') e X 5+(a;), x = 0, 1. (3.6) 

[Hi] Positivity : K{u)' lo) > Q, for almost all {uj,u>') G S-{x) x S+{x), a; = 0, 1. 
(iv) The operator IC{(,,e) is a compact operator from 1/^(5+) onto L^{S-). 

The relation ^3.5|l expresses the conservation of the normal flux of particles at the boundary. 
The reciprocity relation l|3.6> is a macroscopic effect of the time reversibility of elementary particle-surface 
interactions. As a direct consequence of the flux conservation and reciprocity relation, by the change of lu into 
—uj in 13. 5> and the use of 13. 6> . wc get the following 'normalization' identity: 

/ K{uj' ^uj)\uj'Jdiu' = 1, x = 0,l. (3.7) 

The assumption (iv) can be viewed as a regularity hypothesis for the integral kernel K. The adjoint of IC((,,e) is 
obviously given by 

JC*{(t>){x,uj)= [ K{x,uj -^uj')(j}{x,uj')\uj'Jdiu', uieS+ix), x = 0,l, (3.8) 
and is also a compact operator. 
3.1 Flux conservation 

Our first objective is to understand what the effects of the boundary condition H3.H arc on the fluxes of particles. 
For this, we need to define the incoming and outgoing normal fiuxcs associated to distribution function /. 



J-(C,£,i) = / f{^,x = 0,\v\Lu\)\v\\Lu.,\cLj 

J+{i,e,t) = f f{^,x = l,\v\^'\)\v\\u;'Jdu;'. 



li^'es+(x) 

We begin with the following the required identity: 

Lemma 3.2 Assume B('y{f°')) = 0. Then the normal flux is conserved: 

J- = J+. (3.9) 
Proof. The incoming normal fiux is obtained by integrating the left-hand side of the equation 

7"(n=^(7+(n). (3.10) 

Indeed, the coarea formula implies 



J-{e) = [ y-if'^)isiu)V2~e\uj,,\du 

JiueS-{x) 

J+{e) = f -/-{f^){eu;)V2^\cu.,\dco. 



Hypothesis 113.51 and Fubini's Theorem lead to 



J-{e) = [ if K{uj' ^Lj')f{\v\cj'\M\dA\v\\cu.,\du: (3.11) 

JuieS-{x) \Jui'es+{x) j 

f{\v\J\)\v\\J.Ad^' = 3t{e). 

uj'es+(x) 

Finally, it is easy to check that the fiux particles through r± vanished. 

The hvDothcsis l3.1l is crucial for the validity of Theoreni l2.1l From this, the following inequality bears similarities 
with Darrozes-Guiraud inequality in gas dynamics, which we can deduce from Jensen's inequality: 
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Lemma 3.3 Let j+{f) £ L'^{S+) and = X;(7+(/)). Then 

[ |7-(/(.x,a;))|2K,|da.< / \j+{fix,u;))\^u;'Jduj'. (3.12) 

Proof. Using Cauchy-Schwarz inequality and normalization identity 13.711 . we have for x = 0,1: 

|7-(/(x,a;))|2|a;,|d^ = / \lC{-t+{f{x,u;)))\^\u;Miu (3.13) 

S-(x) JS-(x) 



S+{x) 



f K{x,Lu' ^cj)j+{f){x,Lu')\u;',\du;' 

Js+(x) 



{x) 

7+(/(.x,a;))|V.|rf^' 



\uJx\ doj 



Therefore, the operator norm of K, in the space £(L^(iS+), 1/^(5+)) of bounded operator from L'^iS^) to L'^{Sj^) 
is less or equal one: ||A]||£(^2(5^)j < 1. Furthermore, from 13.711 . if ip is constant over 5+, then K.tp is constant 
over 5-. Hence, \\K:\\c(l^i^s+)) = 1- 

We introduce the specular reflection operator J which operates from L^{S+) to according to 

J'ip{u}) = ipiui*), with ujt = (— Wj:, Wy, Wz). Its adjoint J* is defined from (j> e Li^iS-) to 1/^(5+) and is also the 
mirror reflection operator, and we have the equalities J* J = Is_^ = J J* = Is_- Therefore the operator ICJ* 
and its adjoint IC*J operate on I/^(5_) while IC*^ and ICJ* operate on I/^(5+). 

In the derivation of the SHE model, we are interested in the characterization of the equilibria of operator 
/ — ICJ* and I — ]C*J in L'^{S±). By this hypothesis (iv), the operators I — K.J*, I ~ IC*J are Fredholm 
operators. Using Krein-Rutman's Theorem and Fredholm's theory [23, we have the following lemma which can 
be easily adapted from that of |14) . so the proof is omitted: 

Lemma 3.4 (i) The null-spaces N (I — JK.*) and N{I — J*K.) are spanned by the constant functions onS^ and 
Sj^ respectively. 

(a) K.{x,e) IS of norm 1, i.e. \\IC{x,e)(p\\c{ms+)) < Mc(lhs+))- for all ip e L^{S+). 

(m) The range R{I-K.*J) is such that R{I~K.*J) = N{I-J*K.)^. Equivalently, the equation (I-]C*J)f = g 
has a solution f if and only if / g{u})\u}x\dLu = 0. Then, the solution f is unique under the additional constraint 



f{uj)\uj.j,\duj = 0. 



Property (ii) is reminiscent of Darrozcs-Guiraud inequality for boundary conditions in gas dynamics. This 
property expresses that IC{^,e) has a 'good' diffusion behavior. From 13. H , K.ip{ij) for e 5_ appears as a 
convex mean value of p{u)') over 5+. 

Lemma 3.5 The operator K. satisfies 

K.Q+ = Q-K. = JQ+ = JQ-, K.P+ = p-K.. (3.14) 

The proof of this lemma is a straightforward adaptation from Ref. |14) . 
We close this section by giving the following assumption: 

Hypothesis 3.6 There exists ko < 1 such that, for \v\ e R+, ^ e 

ll^-P+ll£(L2(s+),L=(S_)) < fco < 1- 

This assumption comes from the elementary operator theory. Indeed, for every ^ e and |f | £ R^, there exists 
fc(5, |t;|) such that A^P+|| < If] < 1. Obviously, the hvpothesis 13.61 is satisfied in the case of the isotropic 
scattering. 



7 



4 Formal derivation of the macroscopic model 



As a further purpose, we consider a sequence (f")a of solutions to tlie problem (P). Our goal is to study the 
asymptotic behavior of solutions (P) as q — > 0. In this section, we assume that / as q — » in a smooth 

way (in the sense where we will precise later). This approach will formally give rise to a set of diffusion equations 
to be viewed as a SHE system. In fact, this formal derivation will not be completed in this section since an 
'auxiliary problem' - the resolution of which will be postponed to Section El will arise during this study. 

We divide the formal asymptotic into several steps. First, we prove that the limit / is a function of total 
energy. Then we prove that the limit satisfies a continuity equation. An auxiliary problem will be considered 
which allows us to finally derive the current equation. 

The following change of variables will be useful in the remainder of the paper. Since the velocity v satisfies 
the relation. 



\v\ = \/Vx + '"v + '"z, ^^^jii' i = a;,y,z, (4.1) 



we can parameterize the sphere in the direction lj = (ujy.u;,), and ui^ is defined by = cr^l — o;^ — o;^, cr = 
±1. Thus the sphere is given by two local maps {{a,u!y,u:>z), a = ±1}. 
Hence, we have 

Lemma 4.1 Let the change of spherical coordinate be: 

{Vx,Vy,Vz)l i {v,UJx,l^y,(^z) (4.2) 

/ V V 

where uj^ = sgnivx)^ 1 — ujy — uj'^, Uy = y^, ujz = j-^- Then, for sgnivx) = ±1, we have 

df df df uJxOJy df LOxi^z 

■ ^x 



dvx d\v\ du)y \v\ dujz \v\ 

K = (43) 

dVy d\v\ ^ dujy \v\ duo^ \v\ 
K = df 1-iol 

dv^ d\v\ ^ du)^ \v\ du)^ \v\ 



Proof. - For sgn(i;a;) = ±1, we have 



Therefore, for j = x,y,z, we get 



Using 14.41 and the relation 



9/ 


df dv 
dv dvj 


df dujy ^ 
dujy dvj 


df duj^ 
dujy dvj 




d\v\ 












T yjy 




= UJj. 


v\u}j, 


we find 








duJy 






du)z 


UJyUJ^ 


dvy 


\v\ ' 




dVy 


\v\ ■ 


duJy 


UJyUlz 




dujz 




dvz 


\v\ ' 




dvz 


\v\ ■ 


duly 


UlyOJx 




duJz 




dVx 


\v\ ' 




dvx 


\v\ 



(4.4) 



which proves the result. 
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4.1 The limit is a function of the energy 

Taken formally the limit a ^ in ^2. 6(1 - 12.81 shows that / is the solution to equations 

d 



v.Q^'{v><B)-V^jf = 0, (4.5) 
l-{f) = ^(7+(,f)). (4.6) 

We show the 

Lemma 4.2 The solution of \4-5}j - JJ^ are functions of total energy only: 

f{x,v,t) = F{^,^\v\^t). (4.7) 

Proof. The fact that the limit of is a function of the energy only is an easy consequence of Lemma 13.41 
According to Lemma f4. II Eq. 14.51 is equal to 

^{x,u;) + ^^{e.,XLo) = (4.8) 

df 

where — — denotes the derivatives of / with respect to a; G of degree 1, a x b stands for the cross product 

of two vectors a and b and ex x to denotes the tangent vector to in a;-axis. We recall that ^ and v axe mere 
parameters in problem H4.5|l and will be omitted in the remainder of the proof. For f{x,u)) we introduce the 
following change of variables: 

e = -f^ = ^, ujx = aJl-u^l-u^l, a £{-1,1} 

\v\u)x \/2euj V .^g^ 

(.<;* = ijjy cos{6x) + uj, sai{6x), 
uj* = ujy sin(fe) + u>^ cos{9x). 

Since ui = {uiyjCu^), with 14.91 we can write, a;* = {lj*,luI) = 7^+^ (a;), where 

, , /coste — sinfe\ in\ 
K.M^=[sinex cosex) (4-10) 

is a rotation of lo around the x— axis by an angle 6x. Hence the function associated to this transformation is 
written as 

r(a;,a,a;) = /(x,a,7^+,^)(a;)). (4.11) 
Therefore, if / is solution of 12.61 . / satisfies the problem: 

\v\u;^^=0, 7-(,r)=^(7+(.r)). (4.12) 
Integrating the first equation of 1(4.121 with respect to x G [0, 1], yields 

r{x,a,uf) = r (0,(7,6^) = r(l,a,^), (4.13) 
or, in terms of the rotation, TZ^^. ^-^i^^)'- 

fix,a,ui) = r(0,a,7^(,_^)(^^)) = /(l,a,7^+_^.^)(a;)). (4.14) 
This is equivalently written 

7+(,/)=^*7"(/)- (4.15) 
Inserting 14.151 into the second equation of 14.61 leads to 

{I~ICJ*){l-{f)) = (4.16) 
which, by virtue of Lemma [3.41 implies that 7^(/) = -F(||?;p). 
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4.2 The continuity equation 

We arc now concerned with the derivation of the continuity equation I2.12|l . Before, wo introduce the following 
macroscopic quantities. For finite a > 0, wc define the average density F" and current J"(^, e.t) = ( Jy, J") by 

= ^ /"Y r{x,i,\vluj,t)dxdLO, (4.17) 



vs^ 



^'^'^ — 'J2e J J {'Vy(fi){x,V2euj)dxduj, J^{e) = 2e j j unp{x,\/2eio)dxdLo. (4.19) 



r{i,e,t) = ff vr{x,iM,t)dxdLo:=-N{£){vr{e,-)) 

= — I I i^r{x,^,e,uj,t)dxdu;. (4.18) 
a Jo Js^ 

According to reference 1141 . we have, 

Lemma 4.3 Let (fi{x, v) be a function. Then we have 
de 

We prove: 

Lemma 4.4 We temporarily admit that J as a ^ . Then F and J satisfy the continuity equation 12.1^) . 

Proof. Since from Lemma 14.21 / is a function of energy only (in velocity space), it is quite natural to expect 
that be mainly given by a function of e only at order O in a, for instance the macroscopic quantity (/"). We 
expect that there exists a function g°'{^,v,t) such that: 

f''{X,v,t) = F''{^,e.t) +ag"{^,v,t). 

We integrate ^2.6^ with respect to x and on a sphere of constant energy. We get 

v-S/^-E" ■\/.^(F°' + ag"')dxduj = [[ vdcu ■ \'F°' + a\/^ ■ [ [ vg°'dxduj 

-aET ■ [ [ Vyg^dxdiu 



Ana 



(■1 / 



Vj-J^-qS"- / / S/yg"dxduj 



oJs^ 



N{e) 

with V = V{ + E_d/de and for / vdu = 0. Moreover, for any test function ip{e), a straightforward computation 



yields: 



/ i' ( i ^ug°'dxdiuN{e)de = 'in [ i'^-de 

Jc>0 J0JS= Je>0 '^^ 



e>0 JOJS^ Je>0 

which finally proves 

-r— /(ll • - • V^)f"dxdujN{e) = V-r. (4.20) 
Aiva J — ~ 

4.3 The current equation 

Next, we have to find a current equation giving a relation between and the density F°'. We prove: 
Lemma 4.5 The current equation 12. F^) is satisfied. 
Proof. Using Green's formula H5.4^ we get 

d 



:= V. 



dx 



{vx B)- V„/" j xdxdv (4.21) 



v.f'xdv ] U=i - ( / v.rxdv ] |..=o + / / ( v.,^ - (i; X B) • V„ ) xdxdv. 



JR3 
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Now, using 12.81 and the second equation of I2.18|l we get 

\JB.3 J x=l \JS+ JS- 

= \v\'(^J^ 7+(.r)7+(x)k.:|rf^-_^ ^(7+(r))7-(x)k.|rf^j 

= / 7+(r)(7+(x)K|rfu;-^-(7-(x))K.|rfa; = 0. (4.22) 

Similarly, we can prove the relation for (^J^.j v^f^xd^) \^^q- With the first equation of ^2.18|l we get 

/ / .r(vx-^-(vxB)-S/.j}\xdxdv = 2e [ [ uif^xdxdij = aj" . (4.23) 
JoJm-' \ ox ) JoJm'^ 



Therefore, inserting I4.22|l and I|4.2H1 into H4.21|l . wc deduce that 5"* = aJ"" . This justifies the definition of x- 
From irni and lOTt . wc deduce: 



If pi 



'OJR'> JO 

Taking the limit q — > and using 14.711 . wc obtain 



r = -a / dt,rxdxdv- / (v-\/^~\/i:(p-V^) Txdxdv. (4.24) 



J = -[[ (v-S/^- Vjc^ • V„) Fxdxdv, (4.25) 



where J = lim J", F = limF". Taking into account the relation 
we get 

-^^UoL ^' ® ^^^'^'^) (^i " ^i*^!") ^' (^-2*^) 

which leads to equations I2.13|l . Finally, collecting I4.26|l and H4.2()|l . inserting in H4.-5I . multiplying by \/2e and 
taking the limit a — » lead and second to I2.12L I2.17|l . 

It only remains to prove the existence of x which will be achieved in section |^ The relevance of the model 
12. 1211 - 12. 1311 relies on the positivity of the coefficient D. We postpone this task to Section 1^31 

5 Existence of the solutions to the microscopic problem 

In this section, we establish the existence of solutions to the microscopic problem l|2.6> . I|2.8|l . To avoid the 
treatment of initial layers when we pass to the limit a — » 0, we impose a compatibility condition on the initial 
data //: 

Hypothesis 5.1 There exists a smooth function Fj such that fi{x,^,v) = F/(5, |xjp/2) and that fj satisfies: 
fieL'^ie), (v-V^_-V^J-V^fieL\e). 

We define the transport operator by: 

A"f = • - S/^J • V.) / + ^ (^1'.^ - (li X B) . V, j / (5.1) 

on the domain 

DiA"-) = {f e L\e), AVeL^iQ), ^+{f)eL\T+), I'if) = ICj+{f)}. 
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Wc shall denote by A the bare differential operator l|5.1|l when no indication of the domain is needed. Following 
1^ we define the space 

H{A") = {feL'iQ), A'-feL^Q)}. (5.2) 

According to [1], it is well known that the regularity / e H{A") is not suflicient to guarantee the integrability of 
|7+(/)||2(p^) and l7^(/)li2{r_) o^cr the boundary. But if one of these traces is integrable, the other one is also 
integrable. Therefore, we define the following space 

H^{A) = {/ e H'^iA), j-if) e ^^(r^)} = {/ e H^iA), j+if) e L^(r+)}, 

endowed with the graph norm and family of semi-norms, for R> 

=l/ll.,e, +l-4"/ll.<e,' \f\hM = l^^^^^J-^Mf?dT. (5.3) 

Note that hypothesis 15.11 implies that // e D(A°') for all a > 0. For the solvability of the transport equation, we 
require the following hypothesis: 

Hypothesis 5.2 

(i) The first and second derivatives of the potential belong to the Sobolev space W^'°° (in other words, V^(j> is 

bounded and globally lipschitz over R^J. 

(it) The magnetic field B = B{§) e n Wi'°°(R|). 

We give some results concerning the trace operators which are necessary in our setting. Since we will need to 
handle nonlinear functions of / and their traces on the boundary, we have to study for which nonlinear functions 
of / traces can be defined. The problem of existence of a trace is fundamental for the Cauchy problem with 
boundary conditions. This problem was investigate by many authors, such C. Bardos [l]. M. Cessenat in |10| 
studied it and applied it to neutron transport equation, also S. Ukai in (281 for the free transport operator, E. 
Bcals ct al in |S| for the abstract time-dependent linear kinetic equations and recently S. Mischler in I26| for 
Vlasov-Boltzmann equation. 

According to Ref. |S| and hypothesis 15. 21 the following Green's type identity can be easily deduced. 

Lemma 5.3 For f and g in Ho[A") compactly supported with respect to v, we have: 

{A'\f,g)e + {f,A"g)e = - ((7+(/), 7+(ff))r+ - il' {f),l- {g))r^) ■ (5.4) 
a 

From that lemma wc will deduce several estimates for the operator 15. II and their traces on the boundary. 

5.1 Resolution of an approximate kinetic problem 

Applying Leray-Schaudcr's fixed point Theorem to the transport operator A" with domain D{A°') would be 
enough to prove the existence for the kinetic problem. Nevertheless, since the operator IC acts like the identity 
on C*, wc are lacking some estimates on the traces. Following 1141 . we introduce an approximate problem. Let 
?y > be a small approximating parameter; we approach the boundary operator by 

ICr, = P+K. + -^JQ+. (5.5) 
1 + ?7 

We consider the associated iiutial value problem 

adtf:; + Af^ = 0, A;|t=o = F.,„ (5.6) 

E'^ = -V^r, -A^4>"= [ f^dv-C{Q, (5.7) 

l-{fn)=ni+Un)) (5-8) 

with domain D[A'^^) = {/ e H{A"), 7+(/) £ L'^{T+), 7-(/) = ^,,(7+(/))}. We denote by A"^ the transport 
operator A on the domain _D(^^). Thanks to the Lemma f3. 51 we have 

ll^,P+||£{L=(5_),L^{5+)) < 1, V77 > 0. (5.9) 

We can now easily adapt Proposition 4.1 and Lemmas 4.2 and 4.3 in Ref. |14) . Namely we have: 
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Proposition 5.4 (i) For all rj > 0, and for all F,j e D(A'^), there exists a unique function f^ e 
C{[0,T]-D{A^))nC^{[0,T];L'^(e)), that solves fO|) -f5l|). 

(ii) We have also \f^\me) < \Fr,\LHe), \adtf^\LHe) = M"/,"|L^(e) < M"i^,^U2(e). 

( Hi) Let Fj be as in hunothesis \5.1\ There exists a sequence (F,,),,>o such that Fr, £ F>{A'^) and F,j Fj, AF,f — > 
AFj in L^{Q) weak star, as r/ 0. 

5.2 Traces estimates and existence for the kinetic problem 

In this section, wc assume that boundary conditions of tlie form 12.811 are satisfied and prove the existence of 
the traces on the boundary under suitable assumptions. We first estabhsh that the projection P+ of tlic trace 
at the boundary of a function of D{A°') is controlled by the graph norm. This is the key needed to prove the 
convergence of a solution to the problem (P). 
We prove 

Lemma 5.5 We assume that <f) satisfies hu'Dothesis \5.'A If f ^ D(A''^^), then there exists a constant C > such 
that 

\P-l-(f)\\,_^ < \P^1^U)\\,^^ < T^('4'^/'/)e =^ Ca\f\\^ (5.10) 
where ko is such that \\K.P^\\c(l'^{s+),l^{s+)) < '^o- 

Proof. The outline of the proof, which is analogous of that of Lemma 4.4 in is as follows. We apply 
Green's formula 15.41 to f = g times a cutoff function pr{\v\'^/2) = p(\v\'^/2R)) with p e C°°(R+) such that 
< p < 1, p(u) = 1 for u < 1 and p(u) = for u > 2. We obtain, thanks to Lemma r3.4l and hypothesis 13. 61 

2a(A'^pj,,,p„,f)e > (1 - kl) [ |t;,.|P+7+(/)|>Hpdr. (5.11) 

Using the boundedness of 4> leads to 

2a{A''pB.,PRf)e < \A"f\L'.(e) + \f\me) + ^l/lL^(e)- (5.12) 

Finally, taking the limit _R ^ oo in this estimate and recalling 15.1111 allows us to conclude. 

Next, wc give some a priori estimates on the projections C*. Note that, for / e D{A'^), we have 7^(/) = 
/C,j(7+(/)). By Lemma 13.14t and using the form 15. 5> of IC,/, the orthogonal projection of the trace on reads 
Q-'y-(f) = j^^JQ+l+if). Therefore, there exists a function q{f) = q{x,^,\v\), x = 0,1, ^ e \v\ > 0, 
such that 

q = {l + ri)Q-j-{f), on r_, ? = Q+7+(/), on r+. (5.13) 

We deduce 

Lemma 5.6 /// £ D(A'^), then there exists a constant C > such that 

\<l{f)\lR<C{a\f\^^.+R\f\L-He))- (5.14) 
where, for 7i > 0, | • |p jj is defined by \ip\p = \v\<r. h'xWfl'^dT . 

Proof. The key of the proof is the control of Q^^^(f) = Q^7^(/). Notice that, if / G D(A°'), then we can 
straightforwardly verify that sgn{vx)({x)f £ Ho{A°'). We multiply A^f by p ^sgn(vx)({x)f with ( a function 
such that (^{1) = 1, ("(0) = —1 and integrate on F. Thus we apply Green's formula ^5.4^ and take the limit 
R ^ 00 just as in the proof of Lemma [5.51 The remainder of the proof follows now straightforward from 1141 . so 
we omit the detail. 

The existence for the kinetic problem can now be stated. It is given by following proposition, the proof of 
which can be found in Ref. I14| . 

Proposition 5.7 Under hmothesis \H.1\ and \5.i!l there exists a solution f" to problem (P), such that f" G 
L-(0,r;L2(e)), e L-(0,r;L2(e)), P+7+(r) e L°°(0,r;L2(F+))^O+7+(r) e L-(0, T; i^(F+)), for 

all R > 0, and the boundary condition is satisfied in the sense that P^"f^(f°') = SP^l^if")', Q^l^if) = 
(5+7+(/"). Moreover, we have 

r \P+l+{rm\h(r,)dt < C7aV/li=(e)- (5.15) 
Jo 
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6 The rigorous derivation of the SHE model 



In this section, wc prove the convergence of the solutions of the Boltzmann-Poisson equation towards solutioirs 
of the SHE model (with coupled energies). The proof naturally falls into five steps. First, we recall the 
estimate obtaining in the previous section; secondly, we prove the weak convergence of /" and establish the 
continuity equation. Thirdly, wc prove that the current actually converges. Fourthly, we pass to the limit in the 
nonlinear term and we finish the section by investigating some properties of the diffusivity of D. 

6.1 L-^-estimates 

We summarize the estimates deduced from the previous section. In the following estimates, C is a bound for 
different constants depending only on the initial data. 

Lemma 6.1 The following estimates are satisfied by solution f" of the kinetic problem (P) constructed in the 
previous section under hmotheses \3.1\ and \5.<i 

l/"lco(o,T;L2(e)) < \Fi\L-^e), (6.1) 
l^"ric"(o.T;L^(e)) < \A"Fi\L.^e), (6.2) 

/ \P^^^{n\hir^)dt < Ca2|F,|i.(e), (6.3) 
Jo 

I \P-J-{n\Ur )dt < Ca^lF.II (6.4) 
Jo 

f\l{n\l.Rdt < C„\Fi\l^^^^^, (6.5) 
Jo 

where C denotes generic constants independent of a and of the data and q{f") is defined by I5.iy\i> . 

As a straight consequence of H6.1|l . (/")« admits a subsequence (still denoted by (/")) that converges to a 
function /" in L°°{Q,T;L'^{M? x (0,oo)) weak star as a — > 0. From the inequalities 16.311 and 16.411 wc deduce 
that the convergence of the traces P+7+(/") (resp. P-7-(/")) to in ^^(o, T; L2(r_^)) (rcsp. L'^{G,T- L'^{T _))) 
strongly towards zero. Furthermore, using the diagonal extraction process, 16. 5> shows that q{f") converges to a 
function q{x,£^,\v\,t) with a; = 0, 1 in 1/^(0, T; i^(7 x Br)) weak star for any R, where Br is the ball centered at 
and radius R in velocity space by properties of H{div) spaces. In order to give a precise meaning to the limits 
of traces, we notice that the boundedness of A"f" in H^^ (0,T; {&)) by ^2.6^ : this implies (according to 0, 
|18p the boundedness of the sequence {Vxf'^)a\r in H^^O.T; H^/^{j x Br). Thank to this estimate, we obtain 
the convergence in the distributional sense of the traces of /" on F x R'^ to the traces of Finally, the traces 



7±(/'') on r x R| satisfy P-7-(/'') = P+7+(/") = 0, Q--/-{f°) = Q+^+(f°) = q, so that 



f\r = q, (6.6) 
is independent of the angular part of the velocity variable. 

6.2 Weak convergence of /" and the continuity equation 

From Green's formula 15. 4> wc deduce the following weak formulation of Eq. I|2.6> - I|2.8|l : 

Lemma 6.2 Let f" be a solution to the kinetic problem (P) given by Proposition \5.7\ Then, Vt/" £ C^{Q' x 
[0, r],]R^) (i.e. twice continuosly dijferentiable and with compact support znR|x]R+*x[0,r[), such that i>{T,-,-) = 
0, we have 

foL -^^ ("S^ + ■ " - ^-^^^ + ^ (^^-^ " ^ • 

+aj^fiij{0,^,v)de = ^(^J^ £ \v,.\^+{f^h+{i,)~ B*^-{i,)dtdry (6.7) 

Proof. Equation 16.71 is established by writing the weak formulation of the perturbed problem I5.51 - H5.8t and 
letting r] —t in this weak formulation. 

As a consequence, we can prove that /" is a function of the total energy only. 
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Lemma 6.3 The limit function /" is a function of {i^,\v\,t) only, = /"(^j t). 
Proof. Using ^6.7|l with a test function ijj such that 7^(1/') = 0, wc get 



/ / f'^dtde + a [ [ f'^iv-'^S.tp- E°'V^i>)dtde (6.8) 

Jo Je at Jg Jq 



T r , Q 



f" ^'-TT^' - (1^ X 5) • dtde + / fii,\t=o de = 0. 

10 Je \ ox -J Jq 

Hence, when q in 16.811 . using the fact that /" is bounded in L°°{Q,T\ 1/^(6)), we get 

i Iq (""^'^' " ■ ^ii^) '^^'^^ = °- (6-9) 

Eq. 16. 9> supplemented witli 16. 6L proves that /" is a distributional solution of equation 

(^«.:^-(lixB).V„j /'' = 0, f\ = q. (6.10) 

The remainder of the proof is a straightforward adaptation of Lemma [4.21 where wc replace 14. 6t by 7^(/) = ?. 

It remains to prove that Eqs. ^2.12|l - ^2.15|l are satisfied in a weak sense. This is the object of the following 

Lemma 6.4 For any test function e, t) e C^(0' x [0, T]) with support compactly supported m R|xR+x[0,r[, 
we have 

Un^e F"^ + J-(v,- E^) A dtde' 



+ 



/ A■K^/2e Fi^}\t=o de' = 0. (6.11) 



Proof. We apply 15.41 with tp as a. test function. Since tp is even with respect to v, it is easy to check that the 
right-hand side in 16.71 vanishes. On the one hand, 

/ r^dtde^ ( ( f'^ dtde (6.12) 
Jb dt Jo Je dt 

for f"' ^ F in L°°(0, T; 1/^(9)) weak-star and dti' £ L^{'d,T;L'^{Q)). Moreover, on the grounds that 
/e //V'(0,Cil^) = /k2xR+* 47rV^-F/V'|t=o de' an d {vx -§^ ~ {vx B) • V^i/i = 0, the result obviously follows 
from the coarea formula. The continuity equation 12.121 follows from 16.1111 by taking the limit a — > 0. 



6.3 Equation for the current 

Our goal in this section is to prove that the approximate current J_" converges weakly towards a current J that 
satisfies (in a weak sense) the current equation H2.161 (with F the limit of (/")a). We are first concerned with 
the consideration of the most general auxiliary problem of which 12.1811 is a particular case. Given a function 
g{x,uj), find x[x,u)) such that: 

("""^ + X S)V„j x{x,u) = g, {x,u,) e [0, 1] X §2 (6.13) 

7+(x) = ^*(7"(X)) (x,uo)^S+. (6.14) 

Once X is determined, we deduce as a corollary, that there exist functions Xi(2;, tt); ^, e), (i = y,z), solutions of 
problem l|6.13l - H(114|l with the right-hand side g = tUi, unique up to additive functions of ^ and \v\. 

The existence of x is provided by the following lemma which is shown by adapting the proof of Lemma 5.4 in 
|14) in a fairly straightforward way. 
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Lemma 6.5 For all g £ i^([0, 1] x S'^), the problem \(i.l'-!\) - \6-14^ has a solution if and only if g satisfies 



g(x,Lu)dxduj = 0. (6.15) 



Furthermore, if this condition is satisfied, the solution x is unique under the condition 

el 



x{x,Lu)dxdLu = (6.16) 

JoJs^ 

and the set of solutions is the one- dimensional linear manifold {x + F{e), with F{e) arbitrary}. 

Corollary 6.6 The function g = tUi, {i = y, z) satisfies the assumptions of Lemma \(i.,'A 

Thorcforo, the auxiliary function x defined by 12.181 exists and is unique under the constraint 16.1611 . 

We require the following regularity assumptions on Xi- 

Hypothesis 6.7 (i) Xi, {i = y,z), belongs to I/^([0, 1] x for almost every (^, e) £ x R+ and are 
bounded functions on Q away from the set {v.^ = 0}. 

(a) The functions ujiXj{x,uj\^,e) belongs to L^{\f),l] x §^) and I I ujiXjdxduj is a function of {^^e) e 

1R| X R+, uniformly bounded on R| x [0, oo[e and tending to as e ^ 0. 

In order to appreciate the importance of tfiis hypothesis and, in particular, to realize that this hypothesis is 
not empty, we refer the reader to the example of isotropic scattering, where Ref. |14| have computed explicitly 
the diffusion coefficient. 

From Hypothesis 16. 71 (ii) . we deduce that the diffusivity tensor 12.171 is defined and is function of (^,c) £ 
R| X R+. 

We are now in a position to establish the current equation. Actually we shall prove that J" has a finite limit. 

Lemma 6.8 J" is bounded m L^{[0, T]xQ). As a —f 0, J" J in the distributional sense and its limit satisfies 
a weak form of the current equation. More precisely: For any test function tp = {tfiy,i>z) in C^{Q' x [0,r],R^) 



(i.e. twice continuosly differentiable and with compact support in Rj x R+* x [0, r[/. we have 

[ [ ■ i^dtde' ^ [ [ F (s/^- E-^\ {'D'^f)dtd8', (6.17) 



as a 0, where D"^ denotes the transpose o/D. The right-hand side of equation {6.17}) is the weak form of the 
current equation I'A.l'Jd . 

Proof, (i) Since J" is a combination of 7+(J") and 7^(J"), to prove the boundedness of the current in 
L'^([0,T] X Q), it is enough to prove that 7^(J") and 7^ (J") separately are bounded in this space. We prove it 
for 7^(J"), the proof being similar for 7^(J"). First, we note that 

(«7+(,/") = {vP+l+{fn), (6.18) 
by symmetry. Then, using coarea formula and estimate 16.31 we get 



7+(^")li.(lo,Tlxe) = / / -^\{vP+l+{f'^))Ydedt (6.19) 



2 

c_ 

7^ 



dedt 



[ a;P+7+(r)dr 



The result follows straightforwardly. 
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(ii) Now, wc choose (f> as a. test function (j){^,e,t) = \/2£^{^,e,t) ■ x{^,i^',£,t)lp{vx) where lp{vx) is given by 

'^("^■^ = 1 1 \v.A>2p. (^-^O) 

The livpothesis 16. 71 makes it possible to pass to the hinit p ^ in 16.201 . Because of 16. 1311 - 16. 14> we have for 
i e {y,z}: 

Vx^-{vxB)-\/A{V2^xPrXi) = V2^^p^(vx^-{vxB)■\/A{uJ^) (6.21) 



OX -J \ ox 

= -V2eipiXi, in 9 

and 

- lC*-f-{V^i>^x^) = V^Ml+{x^) ~ ic*j-{x-d] = 0, on r. (6.22) 

So tliat, usiiiK <0T|l - <02|l together with coarera formula p. Ill , we get 

/" ( vx-i- -{vxB)- {V2^^^ x) dtde (6.23) 



a ./() Jb V 9x 

K>,|r(7+(%/&-X'A)-^*7"(V2i)-X'0) dtdV 



10 Jr+ 

i-T r i-T 



-II r^a ■ ^jV2edtde = - I I r ■ 'iijdtde', 

« Jo Je - Ja Je - 



Jo Je Jo Je 

which justifies the introduction of the auxiliary function x- Thus, the weak forrmilation ^6.7|l yields: 

/ / riPdtde' = a I I s/TePx ■ ^Htde + a I \/2i//X ' ^\t=ode 
Jo JO' JoJe Je 

f"{v ■ V{ - V^(t> ■ V„)(\/2ex • i^dtde. (6.24) 



Now, by letting a — > in 16.2411 . since x £ -^'^([0, 1] x S^), for almost every (5, e) e R| x R+, the first and second 
terms on the right hand side converge to 0. Then 



lim / / rt]]dtde' =11 f"v-\/i:{V2ex-i')dtde~ I I \/i:(P-\/y){V2ex-^j)dtde 
"^oJoJe' ~ JoJe - ~ ~ JoJe - ~ ~ 

= L1-L2. (6.25) 

On one hand, using 14.71 we have 

Li = I I f ■V^{2eui{x-i')) dt.de = I I FV^-((2ef''^ I / ui{xtp)dxdu?\ dtdO 
Jo Je ~ JoJe ~ V JoJs2 / 

= 11 FV^ ■ (D'^ii) dtde'. (6.26) 
Jo Js' 

On the other hand, using Lemma f4. 191 we get: 

L2 = I I f{Vi(t>V.^){s/2~ex-±) dtde (6.27) 
Jo Je 

= 11 FS/ /:(!)■ (V2e I I S/4V2£x-t)dxdLu^ dtde' 
Jo Je' ~ V Jo Js2 / 

• ^ ("(26)3/2 / / ui{x)dxdiJ) dtde' = I I FS/^-^iB'^f) dtde'; 



lo J0' ~ de \ JoJs2 / Jo J Pi' ~de 

so that, combining H6.261 and 16.27|l to conclude. 
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6.4 Passing to the limit in the nonlinear term 

Notice that the hnear term of H(171 converges because of the weak convergence of f". To pass to the limit in the 
nonhnear term of Eq. I|6.7|l . we need the strong convergence of E". For this task, we use the Poisson equation 
to prove some compactness on the electric fields. 

On the one hand, thanks to the classical regularizing properties of the Laplacian, F°'{£,,e,t) lies in 
L°°{0,T;L'^{&')), implies that, since (j)" solves the Poisson equation, belongs to L'=^{0,f;W^'P{Rl)), 

for all e [l,+oo[. In particular V^cj)" £ Z/°°(0,T; VK^'P(R|)), for some p > 3 and Sobolev's imbedding leads to 
Vj^" e L°°(R|). 

On the other hand, by Poisson equation, we have 

= V^A^\F" - CiO) (6.28) 
The relation 16.2811 together the continuity equation give 

dtE°' = V{ ■ J f" vdxdv. (6.29) 

Hence, as J" is bounded in L^(0, T; Z/^(0')) by Lemma 16.81 we conclude that dtF" is bounded in 
1/^(0, T; VK^^'^(0')) and from 16.2911 the regularizing properties of the Laplacian now yield dt4>" bounded in 
L°°(0,T;W'1'2(R2))^ which yields dtE°' bounded in L°° {G,T- L'^iRf)). Thus, as W^''^^) •--> iLC*?). the 
Aubin-Lions Lemma (see |25| 'l yields that the functional space 

M = {E(i L°°(0, T; Wi'2(R^)), dtE e L^(0, T; i^(R|))}, 

provided with the usual product norm, is compactly embedded in L°°(0, T; I/^^^(R|)) and consequently, as E" 
is bounded in Ml, there exists a subsequence of E" which converges strongly in L°°(0, T; L^^^.(R|)). Its remains 
to prove that E is the solution to the Poisson equation. To do so, we notice that, e L°°{0,T; W^'P(R|)) and 
bounded in this space. Then, up an extraction, we have 

^a^^ L°°(0,T;W^^P{Q')) weak-star. 
Thus, it is easy to notice that (j) is solution of Poisson equation in the sense of distribution. 
Hence, we have 

Lemma 6.9 Under hvvotheses \5. 1\ and \5. M let (/") be a family of a weak solution of the system i2.6V - JF7/\) . 

Then, extracting a subsequence, the family {E°') satisfies 

E"~-*E L-(0,r;LL(K|)) strong, 

as a —t 0. 



6.5 Positivity of the diffusion tensor 

We are now going to discuss some properties of the diffusivity of D. We show that D is a positive definite tensor. 
Therefore l|2.121 - H2.16|l are well-posed. To this end, we prove 

Lemma 6.10 The diffusion tensor D defined by \2.11\j is positive definite; this means that, for all ^ e R^ and 
all eo > 0, there exists C = C(eo) > such that: 

2 2 

(Dy,y) = ^Dyy:,F, >c|F|2 = c^y,2, vy,ceR', Ve>eo- (e.so) 

i,j i=l 

Proof. The proof is inspired from 1141 . Let Y = (i/i, j/2) G R^ such that |y| > and ^{x,u) = X]^=i ViXii^j^)- 
By virtue of ^2.17|l . we get 

(Dy,y) = (26)^/2/ / iTx^y^] iyZ^iy^]dxdLu (6.31) 
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Inserting ^2.18|l into Hd.U^ yields 

(DF, Y) = (2e)3/2 x^y^ (""^^'^ ' ^^^^^^^V^ '^^'^'^ (6.32) 

= (2e)3/2(ylO*<E.,$)5. 
On the other hand, Green's formula leads to 

2{VY,Y) = (2£)'/'(^^ \u,\\^-m^dw- \u,\\^+m^du)j (6.33) 

> (2e„)^/2 (|7-Wli=(s_, - |S*7-Wli=(s,)) > 0. 

On the grounds that \lUx\ is uniformly bounded by positive constant, we deduce that D is definite. Next, assume 
that DF ■ y = 0. Then we deduce that $ does not depend on lj. So that ^''*<3> = lj - Y. Since A°*^ does not 
depend on tj, we have • Y independent of lj for Q in unit disk. This is only possible if F = and D is positive 
definite, and since (^, e) i — > D(5, e) is smooth, we deduce that D is uniformly positive definite. This completes 
the proof. 

Acknowledgements. The author would like to express his gratitude to P. Degond for bringing this problem to 
his attention and T. Goudon for their fruitful comments and suggestions. 
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